\(\int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx\) [49]

Optimal result

Integrand size = 38, antiderivative size = 211 \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\frac {(4 a d+b e+2 b d x) \sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{4 d (a+b x)}+\frac {\left (4 b c d+4 a d e-b e^2\right ) \sqrt {a^2+2 a b x+b^2 x^2} \text {arctanh}\left (\frac {e+2 d x}{2 \sqrt {d} \sqrt {c+e x+d x^2}}\right )}{8 d^{3/2} (a+b x)}-\frac {a \sqrt {c} \sqrt {a^2+2 a b x+b^2 x^2} \text {arctanh}\left (\frac {2 c+e x}{2 \sqrt {c} \sqrt {c+e x+d x^2}}\right )}{a+b x} \]

[Out]

1/8*(4*a*d*e+4*b*c*d-b*e^2)*arctanh(1/2*(2*d*x+e)/d^(1/2)/(d*x^2+e*x+c)^(1/2))*((b*x+a)^2)^(1/2)/d^(3/2)/(b*x+
a)-a*arctanh(1/2*(e*x+2*c)/c^(1/2)/(d*x^2+e*x+c)^(1/2))*c^(1/2)*((b*x+a)^2)^(1/2)/(b*x+a)+1/4*(2*b*d*x+4*a*d+b
*e)*((b*x+a)^2)^(1/2)*(d*x^2+e*x+c)^(1/2)/d/(b*x+a)

Rubi [A] (verified)

Time = 0.11 (sec) , antiderivative size = 211, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.158, Rules used = {1014, 828, 857, 635, 212, 738} \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\frac {\sqrt {a^2+2 a b x+b^2 x^2} \left (4 a d e+4 b c d-b e^2\right ) \text {arctanh}\left (\frac {2 d x+e}{2 \sqrt {d} \sqrt {c+d x^2+e x}}\right )}{8 d^{3/2} (a+b x)}-\frac {a \sqrt {c} \sqrt {a^2+2 a b x+b^2 x^2} \text {arctanh}\left (\frac {2 c+e x}{2 \sqrt {c} \sqrt {c+d x^2+e x}}\right )}{a+b x}+\frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+d x^2+e x} (4 a d+2 b d x+b e)}{4 d (a+b x)} \]

[In]

Int[(Sqrt[a^2 + 2*a*b*x + b^2*x^2]*Sqrt[c + e*x + d*x^2])/x,x]

[Out]

((4*a*d + b*e + 2*b*d*x)*Sqrt[a^2 + 2*a*b*x + b^2*x^2]*Sqrt[c + e*x + d*x^2])/(4*d*(a + b*x)) + ((4*b*c*d + 4*
a*d*e - b*e^2)*Sqrt[a^2 + 2*a*b*x + b^2*x^2]*ArcTanh[(e + 2*d*x)/(2*Sqrt[d]*Sqrt[c + e*x + d*x^2])])/(8*d^(3/2
)*(a + b*x)) - (a*Sqrt[c]*Sqrt[a^2 + 2*a*b*x + b^2*x^2]*ArcTanh[(2*c + e*x)/(2*Sqrt[c]*Sqrt[c + e*x + d*x^2])]
)/(a + b*x)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 635

Int[1/Sqrt[(a_) + (b_.)*(x_) + (c_.)*(x_)^2], x_Symbol] :> Dist[2, Subst[Int[1/(4*c - x^2), x], x, (b + 2*c*x)
/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a, b, c}, x] && NeQ[b^2 - 4*a*c, 0]

Rule 738

Int[1/(((d_.) + (e_.)*(x_))*Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2]), x_Symbol] :> Dist[-2, Subst[Int[1/(4*c*d
^2 - 4*b*d*e + 4*a*e^2 - x^2), x], x, (2*a*e - b*d - (2*c*d - b*e)*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a,
b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[2*c*d - b*e, 0]

Rule 828

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Sim
p[(d + e*x)^(m + 1)*(c*e*f*(m + 2*p + 2) - g*(c*d + 2*c*d*p - b*e*p) + g*c*e*(m + 2*p + 1)*x)*((a + b*x + c*x^
2)^p/(c*e^2*(m + 2*p + 1)*(m + 2*p + 2))), x] - Dist[p/(c*e^2*(m + 2*p + 1)*(m + 2*p + 2)), Int[(d + e*x)^m*(a
+ b*x + c*x^2)^(p - 1)*Simp[c*e*f*(b*d - 2*a*e)*(m + 2*p + 2) + g*(a*e*(b*e - 2*c*d*m + b*e*m) + b*d*(b*e*p -
c*d - 2*c*d*p)) + (c*e*f*(2*c*d - b*e)*(m + 2*p + 2) + g*(b^2*e^2*(p + m + 1) - 2*c^2*d^2*(1 + 2*p) - c*e*(b*
d*(m - 2*p) + 2*a*e*(m + 2*p + 1))))*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, g, m}, x] && NeQ[b^2 - 4*a*c, 0
] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && GtQ[p, 0] && (IntegerQ[p] ||  !RationalQ[m] || (GeQ[m, -1] && LtQ[m, 0])
) &&  !ILtQ[m + 2*p, 0] && (IntegerQ[m] || IntegerQ[p] || IntegersQ[2*m, 2*p])

Rule 857

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :> Dis
t[g/e, Int[(d + e*x)^(m + 1)*(a + b*x + c*x^2)^p, x], x] + Dist[(e*f - d*g)/e, Int[(d + e*x)^m*(a + b*x + c*x^
2)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0]
&&  !IGtQ[m, 0]

Rule 1014

Int[((g_.) + (h_.)*(x_))^(m_.)*((a_) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_)*((d_) + (e_.)*(x_) + (f_.)*(x_)^2)^(q_)
, x_Symbol] :> Dist[(a + b*x + c*x^2)^FracPart[p]/((4*c)^IntPart[p]*(b + 2*c*x)^(2*FracPart[p])), Int[(g + h*x
)^m*(b + 2*c*x)^(2*p)*(d + e*x + f*x^2)^q, x], x] /; FreeQ[{a, b, c, d, e, f, g, h, m, p, q}, x] && EqQ[b^2 -
4*a*c, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\sqrt {a^2+2 a b x+b^2 x^2} \int \frac {\left (2 a b+2 b^2 x\right ) \sqrt {c+e x+d x^2}}{x} \, dx}{2 a b+2 b^2 x} \\ & = \frac {(4 a d+b e+2 b d x) \sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{4 d (a+b x)}-\frac {\sqrt {a^2+2 a b x+b^2 x^2} \int \frac {-8 a b c d-b \left (4 a d e+b \left (4 c d-e^2\right )\right ) x}{x \sqrt {c+e x+d x^2}} \, dx}{4 d \left (2 a b+2 b^2 x\right )} \\ & = \frac {(4 a d+b e+2 b d x) \sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{4 d (a+b x)}+\frac {\left (2 a b c \sqrt {a^2+2 a b x+b^2 x^2}\right ) \int \frac {1}{x \sqrt {c+e x+d x^2}} \, dx}{2 a b+2 b^2 x}+\frac {\left (b \left (4 b c d+4 a d e-b e^2\right ) \sqrt {a^2+2 a b x+b^2 x^2}\right ) \int \frac {1}{\sqrt {c+e x+d x^2}} \, dx}{4 d \left (2 a b+2 b^2 x\right )} \\ & = \frac {(4 a d+b e+2 b d x) \sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{4 d (a+b x)}-\frac {\left (4 a b c \sqrt {a^2+2 a b x+b^2 x^2}\right ) \text {Subst}\left (\int \frac {1}{4 c-x^2} \, dx,x,\frac {2 c+e x}{\sqrt {c+e x+d x^2}}\right )}{2 a b+2 b^2 x}+\frac {\left (b \left (4 b c d+4 a d e-b e^2\right ) \sqrt {a^2+2 a b x+b^2 x^2}\right ) \text {Subst}\left (\int \frac {1}{4 d-x^2} \, dx,x,\frac {e+2 d x}{\sqrt {c+e x+d x^2}}\right )}{2 d \left (2 a b+2 b^2 x\right )} \\ & = \frac {(4 a d+b e+2 b d x) \sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{4 d (a+b x)}+\frac {\left (4 b c d+4 a d e-b e^2\right ) \sqrt {a^2+2 a b x+b^2 x^2} \tanh ^{-1}\left (\frac {e+2 d x}{2 \sqrt {d} \sqrt {c+e x+d x^2}}\right )}{8 d^{3/2} (a+b x)}-\frac {a \sqrt {c} \sqrt {a^2+2 a b x+b^2 x^2} \tanh ^{-1}\left (\frac {2 c+e x}{2 \sqrt {c} \sqrt {c+e x+d x^2}}\right )}{a+b x} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.44 (sec) , antiderivative size = 149, normalized size of antiderivative = 0.71 \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\frac {\sqrt {(a+b x)^2} \left (\left (4 b c d+4 a d e-b e^2\right ) \text {arctanh}\left (\frac {e+2 d x}{2 \sqrt {d} \sqrt {c+x (e+d x)}}\right )+2 \sqrt {d} \left (\sqrt {c+x (e+d x)} (4 a d+b (e+2 d x))+8 a \sqrt {c} d \text {arctanh}\left (\frac {\sqrt {d} x-\sqrt {c+x (e+d x)}}{\sqrt {c}}\right )\right )\right )}{8 d^{3/2} (a+b x)} \]

[In]

Integrate[(Sqrt[a^2 + 2*a*b*x + b^2*x^2]*Sqrt[c + e*x + d*x^2])/x,x]

[Out]

(Sqrt[(a + b*x)^2]*((4*b*c*d + 4*a*d*e - b*e^2)*ArcTanh[(e + 2*d*x)/(2*Sqrt[d]*Sqrt[c + x*(e + d*x)])] + 2*Sqr
t[d]*(Sqrt[c + x*(e + d*x)]*(4*a*d + b*(e + 2*d*x)) + 8*a*Sqrt[c]*d*ArcTanh[(Sqrt[d]*x - Sqrt[c + x*(e + d*x)]
)/Sqrt[c]])))/(8*d^(3/2)*(a + b*x))

Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 3.

Time = 0.70 (sec) , antiderivative size = 215, normalized size of antiderivative = 1.02

method	result	size
default	\(\frac {\operatorname {csgn}\left (b x +a \right ) \left (4 \sqrt {d \,x^{2}+e x +c}\, d^{\frac {5}{2}} b x -8 \sqrt {c}\, d^{\frac {5}{2}} \ln \left (\frac {2 c +e x +2 \sqrt {c}\, \sqrt {d \,x^{2}+e x +c}}{x}\right ) a +8 \sqrt {d \,x^{2}+e x +c}\, d^{\frac {5}{2}} a +2 \sqrt {d \,x^{2}+e x +c}\, d^{\frac {3}{2}} b e +4 d^{2} \ln \left (\frac {2 \sqrt {d \,x^{2}+e x +c}\, \sqrt {d}+2 d x +e}{2 \sqrt {d}}\right ) a e +4 \ln \left (\frac {2 \sqrt {d \,x^{2}+e x +c}\, \sqrt {d}+2 d x +e}{2 \sqrt {d}}\right ) b c \,d^{2}-\ln \left (\frac {2 \sqrt {d \,x^{2}+e x +c}\, \sqrt {d}+2 d x +e}{2 \sqrt {d}}\right ) b d \,e^{2}\right )}{8 d^{\frac {5}{2}}}\)	\(215\)

[In]

int(((b*x+a)^2)^(1/2)*(d*x^2+e*x+c)^(1/2)/x,x,method=_RETURNVERBOSE)

[Out]

1/8*csgn(b*x+a)*(4*(d*x^2+e*x+c)^(1/2)*d^(5/2)*b*x-8*c^(1/2)*d^(5/2)*ln((2*c+e*x+2*c^(1/2)*(d*x^2+e*x+c)^(1/2)
)/x)*a+8*(d*x^2+e*x+c)^(1/2)*d^(5/2)*a+2*(d*x^2+e*x+c)^(1/2)*d^(3/2)*b*e+4*d^2*ln(1/2*(2*(d*x^2+e*x+c)^(1/2)*d
^(1/2)+2*d*x+e)/d^(1/2))*a*e+4*ln(1/2*(2*(d*x^2+e*x+c)^(1/2)*d^(1/2)+2*d*x+e)/d^(1/2))*b*c*d^2-ln(1/2*(2*(d*x^
2+e*x+c)^(1/2)*d^(1/2)+2*d*x+e)/d^(1/2))*b*d*e^2)/d^(5/2)

Fricas [A] (verification not implemented)

none

Time = 0.78 (sec) , antiderivative size = 651, normalized size of antiderivative = 3.09 \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\left [\frac {8 \, a \sqrt {c} d^{2} \log \left (\frac {8 \, c e x + {\left (4 \, c d + e^{2}\right )} x^{2} - 4 \, \sqrt {d x^{2} + e x + c} {\left (e x + 2 \, c\right )} \sqrt {c} + 8 \, c^{2}}{x^{2}}\right ) - {\left (4 \, b c d + 4 \, a d e - b e^{2}\right )} \sqrt {d} \log \left (8 \, d^{2} x^{2} + 8 \, d e x - 4 \, \sqrt {d x^{2} + e x + c} {\left (2 \, d x + e\right )} \sqrt {d} + 4 \, c d + e^{2}\right ) + 4 \, {\left (2 \, b d^{2} x + 4 \, a d^{2} + b d e\right )} \sqrt {d x^{2} + e x + c}}{16 \, d^{2}}, \frac {4 \, a \sqrt {c} d^{2} \log \left (\frac {8 \, c e x + {\left (4 \, c d + e^{2}\right )} x^{2} - 4 \, \sqrt {d x^{2} + e x + c} {\left (e x + 2 \, c\right )} \sqrt {c} + 8 \, c^{2}}{x^{2}}\right ) - {\left (4 \, b c d + 4 \, a d e - b e^{2}\right )} \sqrt {-d} \arctan \left (\frac {\sqrt {d x^{2} + e x + c} {\left (2 \, d x + e\right )} \sqrt {-d}}{2 \, {\left (d^{2} x^{2} + d e x + c d\right )}}\right ) + 2 \, {\left (2 \, b d^{2} x + 4 \, a d^{2} + b d e\right )} \sqrt {d x^{2} + e x + c}}{8 \, d^{2}}, \frac {16 \, a \sqrt {-c} d^{2} \arctan \left (\frac {\sqrt {d x^{2} + e x + c} {\left (e x + 2 \, c\right )} \sqrt {-c}}{2 \, {\left (c d x^{2} + c e x + c^{2}\right )}}\right ) - {\left (4 \, b c d + 4 \, a d e - b e^{2}\right )} \sqrt {d} \log \left (8 \, d^{2} x^{2} + 8 \, d e x - 4 \, \sqrt {d x^{2} + e x + c} {\left (2 \, d x + e\right )} \sqrt {d} + 4 \, c d + e^{2}\right ) + 4 \, {\left (2 \, b d^{2} x + 4 \, a d^{2} + b d e\right )} \sqrt {d x^{2} + e x + c}}{16 \, d^{2}}, \frac {8 \, a \sqrt {-c} d^{2} \arctan \left (\frac {\sqrt {d x^{2} + e x + c} {\left (e x + 2 \, c\right )} \sqrt {-c}}{2 \, {\left (c d x^{2} + c e x + c^{2}\right )}}\right ) - {\left (4 \, b c d + 4 \, a d e - b e^{2}\right )} \sqrt {-d} \arctan \left (\frac {\sqrt {d x^{2} + e x + c} {\left (2 \, d x + e\right )} \sqrt {-d}}{2 \, {\left (d^{2} x^{2} + d e x + c d\right )}}\right ) + 2 \, {\left (2 \, b d^{2} x + 4 \, a d^{2} + b d e\right )} \sqrt {d x^{2} + e x + c}}{8 \, d^{2}}\right ] \]

[In]

integrate(((b*x+a)^2)^(1/2)*(d*x^2+e*x+c)^(1/2)/x,x, algorithm="fricas")

[Out]

[1/16*(8*a*sqrt(c)*d^2*log((8*c*e*x + (4*c*d + e^2)*x^2 - 4*sqrt(d*x^2 + e*x + c)*(e*x + 2*c)*sqrt(c) + 8*c^2)
/x^2) - (4*b*c*d + 4*a*d*e - b*e^2)*sqrt(d)*log(8*d^2*x^2 + 8*d*e*x - 4*sqrt(d*x^2 + e*x + c)*(2*d*x + e)*sqrt
(d) + 4*c*d + e^2) + 4*(2*b*d^2*x + 4*a*d^2 + b*d*e)*sqrt(d*x^2 + e*x + c))/d^2, 1/8*(4*a*sqrt(c)*d^2*log((8*c
*e*x + (4*c*d + e^2)*x^2 - 4*sqrt(d*x^2 + e*x + c)*(e*x + 2*c)*sqrt(c) + 8*c^2)/x^2) - (4*b*c*d + 4*a*d*e - b*
e^2)*sqrt(-d)*arctan(1/2*sqrt(d*x^2 + e*x + c)*(2*d*x + e)*sqrt(-d)/(d^2*x^2 + d*e*x + c*d)) + 2*(2*b*d^2*x +
4*a*d^2 + b*d*e)*sqrt(d*x^2 + e*x + c))/d^2, 1/16*(16*a*sqrt(-c)*d^2*arctan(1/2*sqrt(d*x^2 + e*x + c)*(e*x + 2
*c)*sqrt(-c)/(c*d*x^2 + c*e*x + c^2)) - (4*b*c*d + 4*a*d*e - b*e^2)*sqrt(d)*log(8*d^2*x^2 + 8*d*e*x - 4*sqrt(d
*x^2 + e*x + c)*(2*d*x + e)*sqrt(d) + 4*c*d + e^2) + 4*(2*b*d^2*x + 4*a*d^2 + b*d*e)*sqrt(d*x^2 + e*x + c))/d^
2, 1/8*(8*a*sqrt(-c)*d^2*arctan(1/2*sqrt(d*x^2 + e*x + c)*(e*x + 2*c)*sqrt(-c)/(c*d*x^2 + c*e*x + c^2)) - (4*b
*c*d + 4*a*d*e - b*e^2)*sqrt(-d)*arctan(1/2*sqrt(d*x^2 + e*x + c)*(2*d*x + e)*sqrt(-d)/(d^2*x^2 + d*e*x + c*d)
) + 2*(2*b*d^2*x + 4*a*d^2 + b*d*e)*sqrt(d*x^2 + e*x + c))/d^2]

Sympy [F]

\[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\int \frac {\sqrt {c + d x^{2} + e x} \sqrt {\left (a + b x\right )^{2}}}{x}\, dx \]

[In]

integrate(((b*x+a)**2)**(1/2)*(d*x**2+e*x+c)**(1/2)/x,x)

[Out]

Integral(sqrt(c + d*x**2 + e*x)*sqrt((a + b*x)**2)/x, x)

Maxima [F]

\[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\int { \frac {\sqrt {d x^{2} + e x + c} \sqrt {{\left (b x + a\right )}^{2}}}{x} \,d x } \]

[In]

integrate(((b*x+a)^2)^(1/2)*(d*x^2+e*x+c)^(1/2)/x,x, algorithm="maxima")

[Out]

integrate(sqrt(d*x^2 + e*x + c)*sqrt((b*x + a)^2)/x, x)

Giac [F(-2)]

Exception generated. \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\text {Exception raised: TypeError} \]

[In]

integrate(((b*x+a)^2)^(1/2)*(d*x^2+e*x+c)^(1/2)/x,x, algorithm="giac")

[Out]

Exception raised: TypeError >> an error occurred running a Giac command:INPUT:sage2:=int(sage0,sageVARx):;OUTP
UT:Error: Bad Argument Type

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {a^2+2 a b x+b^2 x^2} \sqrt {c+e x+d x^2}}{x} \, dx=\int \frac {\sqrt {{\left (a+b\,x\right )}^2}\,\sqrt {d\,x^2+e\,x+c}}{x} \,d x \]

[In]

int((((a + b*x)^2)^(1/2)*(c + e*x + d*x^2)^(1/2))/x,x)

[Out]

int((((a + b*x)^2)^(1/2)*(c + e*x + d*x^2)^(1/2))/x, x)